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Abstract 

The low-temperature thermoelectric power of the repulsive-interaction one- 
dimensional Hubbard model is calculated using an asymptotic Bethe ansatz 
for holons and spinons. The competition between the entropy carried by the 
holons and that carried by the backflow of the spinons gives rise to an unusual 
temperature and doping dependence of the thermopower which is qualitatively 

similar to that observed in the normal state of high-T c superconductors. 
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Interacting one- dimensional (ID) electron systems generically exhibit spin-charge separa- 
tion , that is, the Hamiltonian separates at low energies into independent terms describing 
the charge and spin degrees of freedom. To date, the thermopower of interacting ID electron 
systems has only been calculated for a few special cases, where the spin degrees of freedom 
reduce to noninteracting spins 0H, or where the contribution of the spin excitations to 
the thermopower is negligible [|IJ. A more general result for the thermopower of such sys- 
tems is of fundamental interest since the entropy carried by the spin excitations represents 
a qualitatively new type of thermopower, distinct from the familiar contributions of charge 
carrier diffusion, phonon drag, etc. Furthermore, Anderson has argued that the physics of 
the Cu0 2 planes in high-T c superconductors is that of spin-charge separation |4j], and it is 
an interesting question whether their unusual normal state thermopower |5|-|8| can be ex- 
plained on that basis. Since a rigorous treatment of spin-charge separation in 2D systems 
is still lacking, it is clearly of interest to investigate the effects of spin-charge separation on 
the thermopower of ID systems, for which powerful methods such as the Bethe ansatz are 
available. 

In this letter, we calculate explicitly the contributions of both the charge and spin excita- 
tions to the low-temperature thermopower of the repulsive-interaction ID Hubbard model in 
two limiting cases: in the strong-coupling limit, and near the Mott-Hubbard metal-insulator 
transition occuring at half filling. In both of these limits, the charge degrees of freedom 
of the model can be mapped onto weakly interacting spinless fermions fI[,3,|9|,|T0f , while the 
spin excitations can be shown to form an ideal semion gas Jll]] at low temperatures, which 



interacts with the charge degrees of freedom via a backflow condition that ensures that the 
electric current which flows in response to an electric field is really a current of electrons, 
which carry both charge and spin. The competition between the entropy carried by the 
charge excitations and that carried by the backflow of spin excitations leads to a nontriv- 
ial temperature and doping dependence of the thermopower. We comment briefly on the 
possible relevance of these results to an understanding of the normal state thermopower of 
high-T c superconductors and quasi-one- dimensional organic conductors. 



The thermoelectric power S is given in the Kubo formalism by 

s _ 1 ( J o °°{J E (0)J e (t))dt fj\ 
T\j Q °°(J e (0)J e ( t ))dt ej' 

where J e and Je are the electric current and energy current operators in the Heisenberg 
representation, \i is the chemical potential, e is the absolute value of the electron charge, 
and (• • •) denotes the thermal average. In a system with spin-charge separation, the energy 
current can be decomposed into a term associated with charge excitations and a term asso- 
ciated with spin excitations, Je = J% + J%i anc ^ the chemical potential can be written as 
A* = + ^si where \i c = dE /dN + 3F c /dN and /x s = dF s /dN, E being the ground state 
energy, N the number of electrons in the system, and F c and F s the free energies of the 
charge and spin excitations. The thermopower can therefore be expressed as S — S c + S s , 
where S c and S s are defined by Eq. ([I]) with J^ s and /x CjS in place of Je and [i. S c and S s 
can be interpreted as the entropies transported separately by charge and spin excitations 
when an electric current flows. 

We specialize our arguments to the Hubbard model of spin- 1/2 fermions hopping with 
matrix element t between nearest-neighbor sites of a ID lattice with unit lattice constant, 
and subject to a repulsive interaction U when two fermions (of opposite spin) occupy the 
same lattice site. Periodic boundary conditions are imposed on the lattice, which consists of 
L sites, and the system is threaded by a time-dependent magnetic flux (hc/e)<&(t). (In the 
following, we set % — 1.) The complete excitation spectrum of this model was obtained by 
Woynarovich |T2j for the case $(£) = 0; we extend the results of Ref. ]I2] to ^ using 



the arguments of Ref. |JL3[| . In the large-L limit, the energy of the system with iV = L — H 
electrons can be expressed as Wl 



H N s 

E(L -H) = E (L) - ]T e c (k h ) + £ ^(A<x), (2) 

h=l cr=l 



where Eq(L) is the ground state energy at iV = L, obtained in Ref. |Ii 
£c{k) = —2tcosk — At J£°dujJi(uj) cos(u sin k)/[u + uexp(uU/2t)], and e s (A) = 
2t J£°dujJi(u)) cos(cjA) sech{ujU / At) / uj . The momentum [defined modulo(27r)] is given by 

m 



H N s 

P = $(t) - ]T Pc (k h ) + £ P S (A.) + V, (3) 

ft=l er=l 

where V = 7r(^ - N/2 + N s /2 + 1), p c (k) = k + 2 J£°dujJ (uj) sin(usmk)/[u + uexp(uU/2t)], 
and p s (A) = — f£°duJo(uj) sin(o;A) sech(o;C//4t)/a;. The parameters fc^ and A CT are real 
numbers satisfying — 7r < A;^ < 7r, — oo < A CT < oo, and can be interpreted as holes in the 
ground state distributions of pseudomomenta and spin rapidities in the Lieb-Wu equations 



14 1; such excitations are referred to as holons and spinous. (In Eq. ([|), we have omitted 
states with complex pseudomomenta, for which there is a finite energy gap [0.) The holons 
and spinons can not in general be regarded as noninteracting quasi-particles since the kh and 
A CT are not free parameters, but are determined by the asymptotic Bethe ansatz equations 
derived in Ref. 1T31: 



H N a Ms 

Lpdh) = 2-nh + $(£) + J2 Qi (sin fa - smk h >) + J2 ^(si^fa - A CT ) - vr ^ sgn(sin - ReA a ), (4) 

h'=l <r=l a=l 

Lj9 s (A CT ) = -2nJ a - £ 6 2 (A CT - sin**) + £ 0i(A (T - A CT ,) - ^ 2tan- 1 " ° (5) 

fc=l <r'=l a=l W 41 

V 2 tan" 1 = 27 rK a + V 2 tan" 1 " - , (6) 

where M s is the number of down-spin spinons, and the two-body scattering phase shifts are 
given by Q%(x) = 2 / °° dusxn(ujx)/ [u + uexp(uU/2t)}, 02(x) = 2tan _1 [tanh(7rtx/[/)]. The 
charge degrees of freedom are specified by the H quantum numbers Ih [defined modulo(L)], 
which are distinct integers (half- integers) for M even (odd), M being the number of down- 
spin electrons. The spin degrees of freedom are specified by the N s distinct quantum numbers 
J CT , which take values in the range -J max , -J max + 1, • • •, J max , where J max = (N-M -l)/2, 
and by the M s quantum numbers K a , which are integers (half-integers) for N s — M s odd 
(even). The Ih and J a represent holes in the ground state distributions of charge and spin 
quantum numbers in the Lieb-Wu equations, while the K a and the complex parameters X a 
describe the string-like states in the spin sector of the model (their relation to the parameters 
of the Lieb-Wu equations is given in Ref. |12[]). The z-component of the total spin of the 



system is S z = (N — 2M)/2 = (N s — 2M s )/2. For the ground state in zero magnetic field, 
the Ih are consecutive integers (or half-integers) centered about L/2 and N s = N mod 2. 

To evaluate the thermopower ([I]) using the above formalism would be quite difficult in 
general because, while the chemical potential can be extracted from the energy spectrum 
(0), the matrix elements of the current operators would have to be evaluated using the Bethe 
ansatz wavefunctions [14j||, which are quite unwieldy. We therefore consider two limiting 



cases where the weakness of the holon-holon and spinon-spinon interactions allows one to 
construct the matrix elements of the current operators explicitly. Both for H <C L and 
for U ^> t, the holons can be mapped onto weakly interacting spinless fermions [p],p|,^,|lO|, 
allowing the matrix elements of J e and J C E to be evaluated; in addition, for both of these 
cases one can show that in the low-temperature limit (&bT <C v s ) the energy spectrum of 
the spin excitations is 



'tt(1 -5) 



2nJ' 



L 



4L 



(7) 



(7 = 1 \ 

where v s is the spinon velocity in the low-energy limit, calculated in Refs. |T5|JI[], 5 = \1 — 
N/L\, and J' a = J a + sgn(J CT )M s /2. Eq. (|7|) follows upon linearizing the spinon dispersion 
relation in the low-energy (A — > ±oo) limit, replacing the spinon-spinon scattering phase 
shifts by their limiting low-energy form Oi(A cr — KJ) ~ (7r/2)sgn(A (J — A a >), and taking 
tan _1 [4t(A cr — X a )/U] ~ (7r/2)sgn(A cr ). The low-energy spectrum (0) is equivalent to that of 



the "ideal spinon gas" described in Ref. [0: for fixed N s , the degeneracies are equivalent 
to those of a system of spin-1/2 bosons, while the behavior of the Hilbert space as N s is 
varied implies that the spinons are in fact semions, the second term in Eq. (^) representing 
the statistical interaction. The special feature of Haldane's ideal spinon gas is that the 
spinon-spinon interactions are described exactly by mean-field theory, so that spin exchange 
processes between spinons are absent |L1]]. Such processes are implicit in Eqs. @ and (J^j, 
which can be thought of as a nested Bethe- Yang ansatz for the spinons. However, spinon 
spin exchange processes have vanishing amplitude in the low-energy limit, so that the low- 
energy spectrum of the spin excitations (0) is independent of the quantum numbers K a 



which specify the spin wavefunction of the spinons. The excitation spectrum (^) implies a 



where n = N/ L. fi s is not to be confused with the spinon chemical potential, which is zero 
since spinons are thermal excitations. 

Because the spinon-spinon interactions have the mean-field form ([?]) in the low-energy 
limit, they do not contribute to the spinon energy current J E . Consequently, J E commutes 
with the Hamiltonian in the low-energy limit, and has the eigenvalue J E = v 2 s P s (assuming 
equal numbers of right and left movers), where P s is the total momentum of the spinon 
gas. To calculate the contribution of the spinon energy current to the thermopower, we 
make use of the Onsager relation between the thermopower and the Peltier coefficient, and 
define II s = TS S = H E + (j, s /e, where U.% = ( J E ) / (J e ) evaluated for a system subject to an 
infinitesimal electric field £ = —L~ l d<&/dt and with VT = 0. The electric field is equivalent 
to a time-dependent shift of the holon quantum numbers Ih in Eq. (^), which leads to a shift 
in the kh distribution and an electric current J e = (J e ) proportional to S. (The system has a 
finite dc conductivity at T > due to activated Umklapp processes which degrade the electric 
current; however, the ratio of current-current correlation functions in Eq. ([!]) is independent 
of the Umklapp scattering rate in the low-temperature limit.) The shift in the kh distribution 
leads via the spinon-holon scattering phase shifts — @2(A CT — sin kh) in Eq. (||) to a shift in 
the A CT distribution. This backflow of spinons results because the electric field really couples 
to electrons, which carry both charge and spin. Both for H L and for U ^> t, one can 
approximate 02(A CT — sin/^) ~ 2 (A cr ) — sech(27rtA -/f/) (2nt sin kh/U); 62(A cr ) is absorbed 
in a redefinition of the spinon momentum at finite H: p s (A a ) = p s (A a ) + HQ 2 {-h- a ) / 'L. 
Summing over p s (A a ), we obtain 



which we evaluate explicitly below, using the fact that the sum over kh is proportional to 



low-temperature spin entropy per site of nk 2 E 



T/3v s |]TT|, from which it follows that 




n(k B T) 2 dv s 
6v 2 dn ' 



(8) 
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J e , while the sum over A CT is proportional to the excitation energy of the spinon gas, which 



can be evaluated via the correspondence with the model of Ref. (TT . 

We first consider the limit 5 = |1 — n\ <C 1, close to the Mott-Hubbard metal-insulator 
transition. This case has been studied previously via a weak-coupling approximation we 
extend the results of Ref. JT| to arbitrary U, and explicitly verify the assertion of Ref. [Q] that 
the contribution of the spin excitations to the thermopower is negligible in this limit. Eqs. @ 
and (HI) implicitly define an energy band e c (k(p)) for charge excitations, k(p) being the inverse 
of the function p c (k). The holons can be thought of as holes in this energy band, which may 
be approximated near the energy minimum at p = it by e c (k(p)) ~ — (p — it) 2 /2\m*\, 
where /i_ = e c (n) is the T = chemical potential in the limit n — > 1~ [14], and 

1 {1 - 2f™dojJ (u)/[l + exp(ioU/2t)}} 2 



\m 



2t l-2j™dujujJ 1 (uj)/[l + exp(ujU/2t)] 

is the absolute value of the holon effective mass [[HJ . Eq. (|]) implies that the holon momenta 
Pc{kh) differ from those of noninteracting spinless fermions by a term which vanishes as H/L, 
N s /L — > 0; we write ph = p c {kh) = 21x1^/ L + $/L + 5pt, where Sph has contributions from 
holon-holon, holon-spinon, and holon-string scattering. Holon-string scattering merely shifts 
the parity of the holon quantum numbers Ih, and is therefore unimportant in the limit L — > 
oo, while holon-holon and holon-spinon scattering give Sph/(ph — it) — » — 41n2 8/Up'Jji) + 
CT 2 as 5, T — > 0, where C is a [/-dependent constant. The shifts of the holon momenta due 
to the interactions thus have a negligible effect on the charge excitation energies in the limit 
5, T — > 0. Because of the vanishing interactions, the current operators J e and J% commute 
with the Hamiltonian in the limit 5, T — > 0, and Eq. (JJ) can be evaluated straightforwardly 
to give 



k 2 T 



m 



lim o 5 c = sgn(l-n)^^, (10) 



where we have used the electron-hole symmetry of the model about n — 1 flT2| , |10|| . The 
corrections to Eq. ( |IDD due to holon-holon interactions are expected to be (9(<5 _1 ). Evaluating 
Eqs. (H) and ® to obtain /i s and n^, we obtain 



where Iq and I\ are modified Bessel functions. S s is negligible compared to S c in the limit 
<5 — > 0. The low-temperature thermopower of the system thus becomes large and positive 
(hole-like) as the metal-insulator transition is approached from n < 1, and has the opposite 
sign (electron-like) as n — > 1 from above. 

In the limit U ^> t, the full n-dependence of the low-temperature thermopower can be 
calculated: When U/t — > oo, e c (k) — > —2tcosk, p c (k) — > k, and the scattering phase shifts 
in Eq. (f|) are 0(t/U), leading to the well-known mapping of the holons onto noninteracting 
spinless fermions in the strong-coupling limit ||||. The spinon dispersion relation now ex- 
plicitly involves a contribution from the backflow of the holon distribution; Eqs. (@)— © give 
v s = (2irt 2 /U)[l -sin(27m)/27m-81n2(t/7r£f) sin 3 im] + 0(t 4 /U 3 ), which is consistent with 
the results of Refs. []15||[||. The dominant contributions to the low-temperature thermopower 
come from S c and /x s /eT, which combine to give (for n < 1) 

j. k%T / 7r 2 cos 7m \ k\T /cos(27m)/n - sin(27rn)/27r?7, 2 + 241n2 (t/U) sin 2 7mcos7m 

c/^oo ~ et \ 6 sin 2 7m/ 3eJ \ [1 - sin(27rn)/27rn - 8 In 2 (t/icU) sin 3 7ml 2 

(12) 

where J = 4t 2 /U is the antiferromagnetic superexchange constant and 5(2 — n) = —S(n). 
The corrections to Eq. ([12]) are 0(k 2 B T / eU) as T — > 0, and come from holon-holon inter- 
actions, higher order terms in w s , and from Eq. @, which gives H S E /T = TT 2 k B T/6eUn. 
In the physically interesting parameter range 100 > U/t ^> 1, Eq. ([12]) implies that the 
low-temperature thermopower is positive for all n < 1 and negative for all n > 1. The full 
temperature dependence of S c [first term in Eq. (|1^)] can be calculated by the method of 
Ref. (neglecting terms of order ksT/U)] as T increases, S c monotonically approaches the 
high-temperature value sgn(l — n){k B /e) ln[(l — S)/S] 0. The full temperature dependence 
of S s [second term in Eq. (|1^)] is an open problem; however, when k E T ^> J, fj, s is dominated 
by the spin entropy and eW E jk B T 1, so that S s — > sgn(n — l){k B /e) In 2, in agreement 
with the result of Ref. 0. Combining the high-temperature results for S c and S s yields 
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the well-known Heikes formula 0. For small hole dopings 8 -C 1, S is dominated by S c , 
which is large, positive, and a monotonic function of temperature. However, for 8 > 1/3, S 
is dominated by S s , and is negative at high temperatures, but with a positive slope at low 
temperatures, implying the existence of a positive peak in the low-temperature thermopower. 
S has the opposite sign for electron doping. 

The magnetic field dependence of S is also readily obtained in the strong-coupling limit: 
a weak field B applied parallel to the chain does not couple to the charge degrees of freedom 
in the Peierls approximation, but the Zeeman coupling leads for /i£> <C ksT <C J to 

AS(B) = ^^^(^-)\ (13) 



e TTV, 



On \knT 



It is interesting to compare the temperature and doping dependence of the thermopower 
in the ID Hubbard model with that observed in the cuprate materials in which high-T c 
superconductivity occurs, which are widely regarded to be quasi-two-dimensional doped 
Mott insulators The in-plane thermopower of the lightly doped cuprates is generically 
positive for hole doping |j| and negative for electron doping |||7|], with a magnitude which 
increases drastically as the nominal concentration of doped carriers goes to zero, in qual- 
itative agreement with Eq. flTU|). Upon further hole doping, the in-plane thermopower of 
the cuprates universally exhibits a positive peak at low temperatures, then decreases mono- 
tonically, often becoming negative at room temperature in the superconducting samples || ; 
the mirror image behavior, with a negative peak at low temperatures, is exhibited f?j by 
the electron-doped superconductor Nd 2 Cu04_ :c F :c . Similarly, the unusual temperature de- 
pendence of the spinon backflow thermopower in the ID Hubbard model leads for 5 > 1/3 
and U 3> t to a low-temperature peak in the thermopower which is positive for hole doping 
and negative for electron doping. Superconducting samples || of Nd 2 _ x .Ce x Cu04 exhibit a 
peak of the same sign as that observed in the hole-doped cuprates, however, which can not 
be accomodated in a model with electron-hole symmetry, such as the Hubbard model. It 
is of interest, however, that the normal-state thermopower of Nd 2 -xCe x Cu04 goes to zero 
as T — > when a magnetic field is applied to suppress the superconductivity f|, which is 

9 



consistent with the linear low-temperature thermopower we calculate. The smallness of the 
isotropic contribution to the magnetothermopower [Q is qualitatively consistent with our 
result for the ID Hubbard model, Eq. flT3|), which is reduced in magnitude by a factor of 



order ksT/J compared to the high-temperature value |T7|]. The spin contribution to the 
thermopower, which takes the value — (fce/e)ln2 in the high-temperature limit, has been 



previously identified |17]] in a family of quasi-one-dimensional organic conductors which have 
been modeled as Hubbard chains with n = 1/2 and U 3> t; from Eq. (|12|), we would expect 
the low-temperature thermopower of these systems to be positive, and a recent measurement 
|18| on a closely related compound indeed exhibits a positive peak at low temperatures. 
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